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New static regular axially symmetric solutions of SU (2) Yang-Mills-Higgs theory are constructed. 
They are asymptotically fiat and represent gravitating monopole-monopole pairs. The solutions 
form two branches linked to the second Bartnik-McKinnon solution on upper mass branch and to 
the monopole-monopole configuration in flat space on the lower branch, respectively. 
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The nontrivial vacuum structure of SU(2) Yang-Mills-Higgs (YMH) theory allows for the existence of regular 
non-perturbative finite energy solutions, such as spherically symmetric monopoles 0, various axially symmetric 
^vq , multimonopoles 0- 0- EI 0] and monopole-antimonopole systems I n the Bogomornyi-Prasad-Sommerfield 

(BPS) limit of vanishing scalar potential, axially symmetric multimonopole configurations are known analytically Q. 
In these solutions the nodes of the Higgs field are superimposed at a single point. Since in the BPS limit the repulsive 
and attractive forces between monopoles exactly compensate, the BPS monopoles experience no net interaction 0. 
Thus, the BPS configuration with multiple node at the origin can be continuously deformed into a system of individual 
monopoles with unit topological charge (see, e.g., [HI El)- However, as scalar field becomes massive, the fine balance 
of forces between the monopoles is broken and there is only repulsion between non-BPS multimonopole solutions 
t—{ , On the other hand, there are axially symmetric saddlepoint solutions of the YMH theory, where the Higgs field vanishes 
^ ■ at several isolated points on the symmetry axis |fj, |7], la] • Simplest such a solution represent a monopole-antimonopole 
(MA) pair, a magnetic dipole 0,0- 

When gravity is coupled to YMH theory, branches of gravitating solutions arise 0, 0, El ■ The lower branch 
emerges from the fiat space configurations as coupling constant a, entering the Einstein- Yang-Mills-Higgs (EYMH) 
equations, is increases from zero. However, there is a difference between behavior of the gravitating solutions with 
ly-j ' a sin gle o r superimposed nodes of the scalar field [l2l and gravitating monopole-antimonopole chains and vortex 
rings 14]. While in the former case the lower branch ends at a critical value, beyond which the core of the configuration 
would be smaller than the Schwarzschild radius |l2t IT^ ] , for the gravitating monopole-antimonopole chains and vortex 
rings the second branch emerges which is extended back to a = . In this limit the configurations shrink to zero 
size and the scaled solutions approach corresponding solutions of pure Einstein- Yang-Mills (EYM) theory El El El . 

It was mentioned that the additional attraction in the YMH system due to the coupling to gra vity also allows for 
bound monopoles with superimposed zeros of the scalar field which are not present in flat space |l3| . However, one 
can conjecture if the EYMH theory allows for further static axially-symmetric solutions representing, for example, 
separated monopole-monopole (MM) pair. Evidently, beyond BPS limit these solutions do not possess counterparts 
in flat space but gravitational attraction may reinforce the effect of the scalar interaction and a bound state of the 
C$ | gravitating monopoles can exist. . 

In this letter we report the existence of one such new solution representing MM pair. On the upper branch it is 
related to the second Bartnik-McKinnon solution with two zeros [l5|. The properties of the gravitating monopole 
pair are compared with those of MA pair, which on the upper branch is linked to the first Bartnik-McKinnon solution 
with one zero, and with solution for the gravitating charge 2 axially symmetric monopole, linked to the extremal 
Reissner-Nordstrom black hole solution. 

In section II we present the Lagrangian of the EYMH theory the axially symmetric ansatz and the boundary 
conditions. In section III we discuss the properties of the gravitating MM pair. 



II. EINSTEIN- YANG-MILLS-HIGGS MODEL AND AXIALLY SYMMETRIC ANSATZ 

We consider the SU (2) Einstein- Yang-Mills-Higgs theory with action 

S= J 2 1Tr(Vn ~ \ T ^ D ^ D ^)-^-^ 2 -V 2 ) 2 ]v^d 4 x (1) 

Here G and A denote the gravitational and scalar coupling constants, respectively, 7/ is the vacuum expectation value 
of the Higgs field and R is Ricci scalar. 
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In isotropic coordinates the static axially symmetric metric reads 0, fili. Il7| 

ds 2 = —fdt +jdr 2 + —dO 2 + - dip 2 , (2) 

where the metric functions /, m and I are functions of the coordinates r and 8, only. The z-axis {6 = 0, 7r) represents 
the symmetry axis. Regularity on this axis requires m = I there. 
For the gauge and the Higgs field we employ the known ansatz Q 

> K x (n) / (n,m) 

= f -idr + (1 - K 2 )d6\ - nsm9 K 3 ^— + (1 - ^4)^^- d<P , (3) 

$ = 77(*i^ m) +$ 2 r^ m) ) . (4) 

where the su(2) matrices Tr n ' m \ Tg n:,n \ and t^,™' 1 are defined as products of the spatial unit vectors 

-(™,m) = (sin(m#) cos(n(^), sin(m#) sm(rup), cos(m#)) , 
-(n,m) _ ^ cos ^ m g"j CO s(ny?), cos(m0) shx(n(p), — sin(m0)) , 
4 n) = (-sin(nv?),cos(ny),0) , (5) 

with the Pauli matrices r a . 

For m = 2, n = 1 the ansatz corresponds to the one for the MApair solutions 0,00], while for m = 1, n > 1 it 
corresponds to the ansatz for axially symmetric multimonopoles |.l la Il3|. In particular, for m = 1, n = 2 we have 
non-BPS extension of the charge 2 monopole solution 

The four gauge field functions Ki and two Higgs field functions $^ depend on the coordinates r and 9, only. To 
construct regular solutions we have to fix the gauge condition rd r K\ — deK 2 — 5]. Further, we introduce the 
dimcnsionless coordinate x = err] and rescale the Higgs field as $ — > Q/r). Then the dimensionless coupling a, 
a 2 = AttGij 2 enters the equations. 

To obtain asymptotically flat solutions which are regular and corresponds to the gravitating MM pair, we need to 
impose the boundary conditions. Regularity of the solutions at the origin (r = 0) requires for the metric functions the 
boundary conditions d r f(r,6)\ r= o — d r m(r,8)\ r = — d r l(r,9)\ r =Q = , whereas the gauge field functions Ki satisfy 
Ki(Q, 6) = K 3 (0, 6) = 0, K 2 (0, 9) = K 4 (0, 6) = 1 , and the Higgs field functions satisfy 

sin<9 $i(0, 6») +cos(? $2(0, 6>) = , (6) 



9 r [cos0 $i(r,0)-sin0 $ 2 (r, 0)]\ r=Q = . (7) 

These conditions are the same both for MA-pair and MM-pair. For the charge 2 monopole both $1 and $2 must 
vanish at the origin. 

The boundary conditions at infinity shall provide correct asymptotic behavior of the fields depending on the con- 
figuration. Evidently, asymptotical flatness requires / — ► 1, m — > 1, I — > 1 for any solution. Also the Higgs field 
at infinity have to approach the 'hedgehog' shape, i.e., <&i — ► 1, $2 — ► . But the boundary conditions on the 
gauge functions can be different. 

To construct MA-pair, which is a deformation of the topologically trivial sector, the gauge field at infinity required 
to tend to a pure gauge — ► id^UW , where U = exp{-i6»4™ ) } H 0. In terms of the functions Ki these 
boundary conditions read: 

K x — >0, K 2 — K 3 — > 2sin6», K t — >l-2cos0, (8) 

whereas for the charge 2 monopole we required 

K x — » , K 2 — > , K 3 — -4 , K 4 — » . (9) 

To construct monopole- monopole pair configuration on the same axially symmetric ansatz 10,®, let us note that 
the multimonopoles can be nicely described in terms of the effective electromagnetic quantities |8|, [18| , like magnetic 
charge 

9 = h I \ Tl ( F *i Dk *) £l ^ r (10) 
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Figure 1: The gauge field function K4 (left) and the metric function of the lower branch monopole-monopole pair solution are 
shown at A = 0.5, a = 0.05 



and the dimensionless magnetic dipole moment /1, which is associated with asymptotic behavior of the function K3 
as K3 — ► (1 — cos 6)/ sin 8 + /is'mO/r For a MA pair, for example, the integrated magnetic charge is zero, however 
the charge density distribution g(x) = ^Tr {FijDk^>) £ijk is not trivial, it has a maximum associated with node of the 
Higgs field on positive z-axis and symmetrically located minimum associated with the node on negative z-axis. The 
MA pair has non-vanishing magnetic dipole moment which can be relatively good evaluated by consideration of the 
magnetic charges as point charges located at the nodes 0, ^| . On the other hand, the charge 2 monopole posses zero 
dipole moment since both nodes coincide. 

Therefore, we can conjecture that the MM pair possesses zero dipole moment. This condition implies that the 
boundary conditions on the functions, which enter the component A v of the gauge potential, have to be modified and 
we impose 

Ki — ► , K 2 — ► -1 , r 2 d r K 3 — ► , r 2 d r K 4 — ► , (11) 

We find that this modification yields new branch of gravitating MM solutions. 

Finally, regularity on the z-axis requires dgf — d$m = dgl — 0, whereas the matter field functions satisfy K± = 
K3 = $2=0, dgK2 — dgK^ — dg&i — for all these configurations. 



III. NUMERICAL RESULTS 



Subject to the above boundary conditions (|llfl . we solve the system of 9 coupled non- linear partial differential 
equation numerically in compact radial coordinate x = r/(l + r) € [0 : 1]. The numerical calculations are based on 
the Newton-Raphson iterative procedure [19|. 

For a small non-vanishing values of A and a we obtain the new static solution with two nodes of the Higgs field 
on the z-axis, which smoothly evolves as coupling constants begin to vary. Furthermore, there is a limiting branch of 
gravitating monopoles as A = 0. These solutions are quite different from the known MA configurations (see Fig^), in 
particular, the metric function / does not possess a minimum at the origin. 

To confirm that these solutions can be interpreted as the monopole-monopole pair, let us consider the magnetic 
charge density defined by integrand in (|10fl . For the configuration which we are analysing, the charge density remains 
positive everywhere. It has a shape of two tori those maxima form two rings in planes parallel to the a;?/-plane, 
intersecting the symmetry axis close to the nodes of the scalar field (see Fig [3 left) . The energy density of the 
configuration possesses two maxima on the z-axis associated with nodes of the scalar field (Fig [21 right). 

We can evaluate the charge by straightforward substitution of the numerical solutions into the definitions above, 
the calculation gives g = 2. Thus, the solution represent gravitationally bounded pair of monopoles. In the limiting 
case a — > the solution approaches flat space limit where the bound state ceases to exist and monopoles set to be 
free. Separation between the monopoles in this limit depends on the value of the Higgs self-coupling A, it is minimal 
for BPS monopoles: d m i„ = 6.53 in scaled units. As a increases the monopoles approach the flat space limit on a 
larger separation, e.g., d m i n — 7.31 for A — 0.1. 
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MM configuration: charge density at X=0.1, (1=0.03 MM configuration: energy density at X=0.1, a=0.03 




Figure 2: The charge density (left) and the energy density (right) distributions of the monopole-monopole pair solution are 
shown as functions of the coordinates z, p = rsinf? at A = 0.1, a = 0.03. 




Figure 3: The nodes of the Higgs field are shown as functions of the coupling constant a for the MM pair and the MA pair 
solutions at A = (left). Also the scaled masses of the MM pair, the MA pair and the charge 2 monopole solution are shown 
as functions of the coupling constant a at A = 0. 

This branch of gravitating monopoles extends up to a maximal value a cr , where it merges with a second branch 
similar to the case of the gravitating MA pair ^4|. As scalar coupling increases the a cr decreases, e.g., for A = we 
have a cr — 0.838 while for A = 0.1 a cr = 0.616 and for A = 0.5 a cr — 0.562. The second branch extends back to 
a = 0, as seen in Fig EJ In this limit the mass diverges and the configuration shrink to zero size. However, rescaling 
the coordinate x — ► ax and the scalar field leads to a limiting solution with finite size and finite scaled 

mass Jjj . The scaled mass of the MM pair and the nodes of the Higgs field as functions of the coupling constant a 
are exhibited in Fig [21 Note that there is a principal difference between the metric function / of the gravitating MM 
pair and that of the MA pair: the minimum of the former function is clearly associated with position of nodes of the 
scalar field while the minimum of the latter function remains at the origin Q ■ 

Considering the limit a — > on the upper branch in scaled coordinates, we observe that the solution may be thought 
of as composed of a scaled charge 2 monopole solution in the inner region and the second Bartnik-McKinnon solution 
with two nodes of the gauge field function in the outer region, as seen in Fig0] The inner region shrinks to zero size 
as a — ► 0. 

Comparing the solutions with the gravitating charge 2 monopole [l3| we observe that on the lower mass branch, the 
mass of the monopole-monopole pair for the same value of A = is lower (see Fig El- Since both configurations are in 
the same topological sector with charge 2, one may expect that gravitating 2-monopole is unstable. Hovewer, these 
solutions evolve differently, as a increases, the branch of gravitating 2-monopole approaches the extremal Reissner- 
Nordstrom black hole with magnetic charge 2 ^2] while both gravitating MM pair and gravitating MA pair evolve 
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MM configuration: f at X=0, u=0.D6 MM configuration: K 2 at X=0, a=0.06 
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Figure 4: The metric function / (left) and the gauge field function K2 (right) of the MM pair solution are shown in scaled 
coordinates at a — 0.06 and A = 0. 
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Figure 5: The metric function / (left) and the gauge field function K2 (right) of the upper branch MM pair solution, the upper 
branch MA pair solution and the charge 2 monopole are shown at A = and a = 0.05 (MM and MA configuration) and at 
a — 1.45 (charge 2 monopole) 



back on the upper branch and they are linked to the Bartnik-McKinnon solutions of pure EYM theory. 

In is instructive to compare the limiting behavior of these three different configurations, all possessing two nodes 
of the Higgs field. In Fig [S] we exhibit the corersponding gauge field functions K 2 and the metric functions / for the 
upper branch MM pair solution, the upper branch MA pair solution and charge 2 monopole. 

Concluding, we have found new static axially symmetric solutions of SU(2) EYMH theory, which represent grav- 
itating monopole-monopole pairs. This results holds both for zero and for finite Higgs self-coupling |20|. We expect 
that there are also gravitating multiply magnetically charged solutions which are counterparts of the vortex rings 
discussed in p| . Since we observe that for the MM pair the minimum of the metric function / coincides with location 
of the monopoles, some interesting features of the black hole solutions related to these configuration may be observed. 
Study of these and other solutions of the new class, only simplest of those was discussed in this note, will be presented 
elsewhere |2fjj. 
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